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A  linear  theory  for  the  transverse  motion  of  a  cable  being  towed  without 
purposeful  maneuvers  was  developed. 

A  cable  experiment  was  then  performed  at  sea  to  verify  the  linear  theory 
and  to  establish  parameters  that  would  apply  operationally  in  the  sea  environ¬ 
ment.  Specifically,  the  linear  transfer  function  relating  the  transverse 
horizontal  motion  of  the  upstream  and  downstream  ends  of  a  long  steel  cable 
was  calculated  and  measured. 

Finally,  the  theoretical  results  were  compared  with  the  experimental 
results.  The  agreement  between  the  results  was  reasonable,  but  not  perfect. 

A  more  exact  treatment  of  the  three-dimensional  problem,  and/or  incorporation 
of  the  usual  nonlinear  hydrodynamic  forces,  would  presumably  improve  the 
model.  However,  the  major  effort  required  would  result  in  only  a  minor 
improvement . 

The  results  of  this  study  are  an  improved  understanding  of  the  transverse 
dynamics  of  long,  flexible  cylinders  in  axial  flow,  and  an  analytically  simple 
(relative  to  numerical  methods)  algorithm  for  predicting  the  horizontal  motion 
of  a  towed  vehicle  caused  by  the  random  crosstrack  meandering  of  the  towing 
vessel. 
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INTRODUCTION 


GENERAL 

Ocean  engineers  frequently  conduct  studies  in  which  a  surface  vessel 
tows  an  oceanographic,  optical,  or  acoustic  instrumentation  package  at  the 
end  of  a  long  cable.  It  is  important,  in  such  a  case,  to  know  the  position 
of  the  instrumentation  package  while  the  vessel  attempts  to  maintain  a 
constant  course.  The  degree  of  precision  needed  when  predicting  this 
position  depends  on  the  task  itself.  This  precision  can  vary  from  tens  of 
meters  in  oceanographic  applications  to  tens  of  centimeters  for  some 
acoustic  applications.  This  report  documents  the  findings  of  a  study 
undertaken  to  develop  procedures  for  predicting  the  position  of  the 
instrumentation  package  towed  behind  the  cable.  This  study  should  be  viewed 
as  an  effort  to  develop  an  analytical  basis  for  an  optimum  linear  shape 
predictor.  Such  an  estimator  requires  additional  information,  such  as 
optimal  parameter  values,  instrumentation  error  characteristics,  etc. 

In  the  past,  approximate  numerical  schemes  have  been  used  to  simulate 
the  dynamic  conditions  of  interest  to  this  study. ^  Numerical  techniques 
require  numerical  experiments  to  obtain  insight  into  the  physical  processes 
involved.  Extensive  experimentation  is  typically  not  possible  because 
numerical  studies  require  significant  computer  capability  and  are  tedious 
and  costly.  Mudie  and  Ivers  concluded  that  sinusoidal  deviations  of  the 
towing  ship  from  a  straight  towing  track  resulted  in  delayed  and  reduced 
excursions  of  the  towed  package. ^  This  result  will  be  commented  on 
extensively  in  the  present  work,  in  which  a  different  approach  was  taken; 
explicit  solutions  of  partial  differential  equations  describing  the  cable 
dynamics  were  sought.  More  specifically,  the  cable  was  considered  to  be  a 
long  flexible  cylinder  in  axial  flow  in  a  two-dimensional  horizontal  space. 
The  consequences  of  the  cable  being  negatively  buoyant,  and  thus  occupying 
depth,  the  third  dimension,  were  handled  in  an  ad  hoc  manner  by  modifying 
the  parameters  to  account  for  the  horizontal-vertical  interaction  of  forces. 

This  study  comprised  an  analytical  phase  and  an  experimental  phase  with 
a  resultant  data  comparison.  In  the  Analytical  Program  section  of  this 
report,  the  momentum  equation  developed  by  Paidoussis^  and  the  particular 
solution  suggested  by  Ortloff  and  Ives^  are  discussed.  Various  aspects  of 
the  Ortloff  and  Ives  solution  are  then  discussed  in  detail;  the  solutions  to 
the  particular  problem  of  sinusoidal  forcing  of  the  upstream  boundary  of  the 
cable  are  discussed  first,  and  then  the  complementary  solution  is 
developed.  The  complementary  solution  leads  to  further  discussion  regarding 
the  stability  of  the  system.  The  section  concludes  by  developing  a  method 
for  modifying  the  results  to  incorporate  the  effects  of  the  non-neutrally 
buoyant  cable.  In  the  Experimental  Program  section,  a  cable  experiment 
performed  at  sea  is  described;  the  experiment  was  designed  to  verify  the 
theory  and  to  establish  parameters  that  apply  operationally  in  the  sea 
environment.  A  theoretical  experimental  comparison  discussion  follows. 
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ANALYTICAL  APPROACH 

Our  study  began  with  an  examination  of  the  behavior  of  long,  neutrally 
buoyant  cylinders  under  tow  in  water;  the  study  was  limited  to  the  case  in 
which  the  cable  assumes  relatively  shallow  angles  with  respect  to  the  flow, 
which  is  typical  when  the  tow  ship  attempts  to  steer  a  straight  course.  We 
were  primarily  concerned  with  the  random  lateral  motion  about  a  straight 
towing  path.  The  aim  of  the  analytical  portion  of  the  study  was  to  solve  a 
momentum  equation  in  order  to  predict  the  position  of  the  instrumentation 
package. 


Paidoussis'  Equation 

The  basic  approach  used  the  solution  suggested  by  Ortloff  and  Ives^ 
to  solve  the  momentum  equation  developed  by  Paidoussis. 2 

We  are  interested  in  a  momentum  equation  that  describes  the  dynamics  of 
a  long,  thin,  flexible  body  in  fluid  flow.  The  equation  is,  necessarily,  a 
nonlinear  partial  differential  equation  as  a  consequence  of  the  nonlinear 
hydrodynamic  forces.^  Solving  this  nonlinear  equation  would  be  quite 
difficult,  and  no  one,  to  date,  has  attempted  it. 

Paidoussis,  however,  suggested  that  the  equation  could  be  linearized 
when  the  cable  assumed  small  angles  with  respect  to  the  flow*.  This  change 
from  the  nonlinear  to  the  linear  equation  greatly  facilitated  the 
development  of  a  solution.  Paidoussis  based  his  linearization  theory  on  a 
modification  of  the  model  for  hydrodynamic  forces  suggested  by  G.  I. 

Taylor^  in  his  work  on  the  swimming  of  eels  and  worms.  However,  the  basis 
for  the  Paidoussis  linearization  is  questionable  because  of  two  factors: 

(1)  this  modification  was  not  developed  for  the  smooth  cylinders  of  interest 
to  Paidoussis  (and  to  this  study),  and  (2)  the  form  of  the  hydrodynamic 
force  at  shallow  angles  that  allowed  for  the  linearization  had  not  been 
experimentally  verified.  Recently,  Rispin^  performed  experiments  that 
verified  the  "small  angle"  form  used  by  Paidoussis  in  his  linearization  and 
defined  the  limits  over  which  such  a  linearization  applies.  Specifically, 
it  appears  that  the  linear  momentum  equation  derived  by  Paidoussis  applies 
if  the  cable-fluid  angles  do  not  exceed  2.5  to  3  degrees. 

The  purpose  of  the  solution,  as  worked  out  by  Paidoussis,  was  not  to 
determine  shape  explicitly,  but  rather  to  predict  the  instabilities  that 
might  exist  in  such  a  system.  Unfortunately,  the  complexity  of  even  the 
linearized  expression  prevents  easy  interpretation  or  understanding  of  the 
physical  process.  Numerical  analysis  of  the  power  series  solutions  was 
performed  by  Paidoussis  to  identify  stability  regions  and  was  later  extended 
by  Pao.^ 


♦The  concern  here  is  for  the  dynamically  changing  cable  angle  in  the 
horizontal  plane.  The  static  vertical  depression  angle  resulting  from  the 
weight  of  the  cable  will  be  addressed  later. 
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Pao's  Contribution 


Pao's^  stability  calculations  are  also  particularly  pertinent  to  this 
current  work.  His  theoretical  and  experimental  work  showed  that  there  are 
no  hydrodynamic  instabilities  in  long  cylinders  because  long  cylinders  have 
tensions  high  enough  to  filter  out  the  modes  of  instability,  except  near  the 
free  end  (within  20  diameters)  where  the  tensions  again  become  quite  small. 
Pao's  results  agree  with  the  tow  tank  observations  at  the  David  Taylor  Naval 
Research  and  Development  Center^,  which  also  showed  that  no  instabilities 
occur  along  towed  flexible  cylinders.  The  conclusion  directed  this  effort 
to  study  motions  in  towed  cylinders  that  result  from  external  forces  on  the 
system,  not  spontaneous  instabilities.  Thus,  we  should  expect  the  cable  to 
respond  to  forces  applied  by  the  tow  vessel  and/or  oceanographic  current 
shears . 


Previous  Attempts  to  Solve  Paidoussis*  Equation 

Although  no  one  has  attempted  to  solve  the  full  nonlinear  momentum 
equation,  several  modest  attempts  have  been  made  to  solve  the  Paidoussis 
equation  with  the  intention  of  studying  cylinder  shape.  Pao  and  Tran^ 
used  numerical  integration  schemes  to  solve  the  Paidoussis  equation  with  a 
sinusoidal  forcing  at  the  forward  end.  Only  numerical  examples  were  worked 
out  and  no  general  conclusions  about  the  system  resulted.  Gardner  and 
LindemannlO  considered  the  solution  for  the  special  case  when  the  ratio  of 
the  forcing  wavelength  to  cylinder  diameter  is  snail  (less  than  50).  They 
showed  that  the  motions  caused  in  this  case  were  the  same  as  in  a  wave 
system  traveling  in  a  "beam";  for  high  frequencies  the  flexural  rigidity  of 
the  cable  dominates  and  thus  causes  the  cylinder  to  behave  like  a  beam. 

They  also  showed  that  the  fluid  forces  were  small  relative  to  the  flexural 
rigidity  of  the  cylinder  so  that  the  actual  response  of  the  system  was 
similar  to  that  of  a  beam  in  a  vacuum.  However,  Gardner's  and  Lindemann's 
solution  is  of  little  interest  to  our  study  because  the  tow  vessel  will 
excite  the  cable  with  longer  wavelengths  as  it  meanders  back  and  forth 
randomly  across  its  average  course. 

An  explicit  solution  to  the  Paidoussis  equation,  modified  by  the 
omission  of  the  flexural  stiffness  term,  was  found  by  Ortloff  and  Ives. 3 
This  solution  applies  in  the  parameter  region  opposite  to  that  addressed  by 
Gardner  and  Lindemann;  the  Ortloff  and  Ives  solution  assumes  that  the 
forcing  wavelength  is  much  larger  than  50  times  the  cylinder  diameter,  thus 
making  the  cylinder  respond  as  though  it  were  infinitely  flexible.  This 
modification  is  particularly  applicable  to  this  study,  in  which  we  are 
interested  in  the  low  frequency  cross  track  meandering  of  a  ship  towing  a 
cable  where  forcing  wavelengths  typically  exceed  1000  meters.  The  principal 
result  of  the  Ortloff  and  Ives  work  was  an  explicit  closed  form  solution  for 
the  cylinder's  space-time  shape. 

Kennedy! 1  developed  approximate  solutions  to  the  Paidoussis  equation 
for  low  and  high  frequency  excitation  of  the  system  by  the  tow  vessel.  He 
determined  that  the  cable  responds  in  three  distinct  manners  over  a  range  of 
excitation  frequencies.  At  low  frequencies,  the  cable  is  dominated  by 
hydrodynamic  forces  that  cause  deformations  put  into  the  cable  by  the 
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towing  vessel  to  propagate  downstream  relatively  unattenuated  and  with  a 
propagation  speed  near  the  tow  speed.  As  the  frequency  increases,  the 
inertial  forces  and  tension  stiffness  of  the  cable  become  important  and  wave 
systems  typical  of  string  waves  are  produced.  These  waves  have  a  phase 
velocity  that  diminishes  as  the  wave  proceeds  toward  the  free  end  as  a 
consequence  of  the  diminishing  tension.  Finally,  at  high  frequencies,  the 
response  of  the  cable  is  dominated  by  the  flexural  rigidity  of  the  cable  and 
the  "beam  in  a  vacuum”  response  noted  by  Gardner  and  Lindemann  is 
reproduced.  This  present  study  goes  beyond  the  approximate  solutions  of 
Kennedy's  work^  to  develop  an  analytically  more  exact  solution. 


ANALYTICAL  PROGRAM 


GOVERNING  EQUATION 

The  basic  partial  differential  equation  (PDE)  applied  in  this  study  was 
derived  and  discussed  in  several  papers. ^ » 3 , 1 1  equation  describes  the 

dynamics  of  the  physical  configuration  (see  figure  1)  that  consists  of  a 
flexible  cylindrical  body  of  circular  cross  section,  immersed  in  an 
incompressible  fluid  of  uniform  density  which  is  flowing*  with  uniform 
velocity  (U)  parallel  to  the  x  axis.^  The  equation  is  expressed  as: 


(M+m)  i_I  ♦  MU2  -Lz  +  2MU  fj§£  -  |j 


_  .2  2c 
+  2c!  MU  -2-Z  +  — 


*  s  [2c<  Ir  u2a-*)  H] 

H  °(&  *  °  *)■  °  - 


where 


y(x,t) 


*  crosstrack  position  of  the  cylinder  measured  at  x, 

*  along  track  independent  space  variable  increasing  downstream, 

■  time, 

*  effective  mass  of  fluid  "pushed"  by  the  cylinder  per  unit 
length  of  cylinder, 

*  mass  of  the  cylinder  per  unit  length, 

■  free  stream  velocity, 

*  cylinder  length, 

■  cylinder  diameter, 

*  tangential  (longitudinal)  drag  coefficient, 

*  drag  coefficient  added  to  account  for  a  drogue  at  the  end  of 
the  cable,  and 

■  coefficient  of  the  linearized  normal  drag  expression. 


*0r  equivalently,  being  towed  with  uniform  velocity  U. 
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Figure  1.  Towing  System  Coordinates 


Equation  (1)  is  the  Paidoussis^  equation  with  the  following 
modifications  and  clarifications: 


1.  The  fourth-order  spatial  derivative  describing  the  transverse  force 
due  to  the  cable's  flexible  rigidity  was  neglected.  As  a  rule  of  thumb, 
this  assumption  is  appropriate  for  a  cylinder  with  a  diameter  (dc)  less 
than  2  percent  of  the  deflection  wavelength. 

2.  The  form  of  the  normal  hydrodynamic  force  (sixth  term)  is  a  result 
of  the  Paidoussis  linearization  discussed  in  the  introduction.  Laboratory 
measurements  by  Rispin^  indicate  that  this  form  is  appropriate  for 

cable/ flaw  angles  of  less  than  2.5  to  3  degrees. 

3.  The  form  of  the  hydrodynamic  drag  coefficients  (ct,  c^,  and 
cn)  was  changed  to  be  the  same  as  that  typically  found  in  modern 
literature  and  defined  by  the  normal  and  tangential  drag  forces  expressed  by: 


Ft  ■  f  dc»Vt  *nd  Ft 


§•  d  c  V2 
2  c  n  n 


(lb) 


where  vc  and  vn  are  the  components  of  U  tangential  and  normal  to  the 
cylinder  respectively;  p  is  the  fluid  density.  The  relation  between  these 
drag  coefficients  and  those  used  by  Paidoussis,  Ortloff,  and  Ives  is 


c 


N 


4 

—  c  and 
IT  n 


*  4c. 


where  cjj  and  ct  are  the  Paidoussis  coefficients. 


(lc) 
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The  following  nond imens iona 1  variables  were  introduced  into  equation  (1): 

T  3  u  u  ’  =  z  *  n  a  a  ’  3  •  e  a  b  »  and  '*> a  2*f  | .  <2> 

c 


where  f  is  a  frequency  variable  and  A  is  an  arbitrary  length  which  can  take 
on  any  value  since  equation  (1)  is  linear  in  y.  (Later,  A  will  be 
identified  with  the  magnitude  of  the  forced  vibration.)  The  result  is 


a2n 

3t2 


+  2S 


a2n 


2 

,  .  _v  9  n  9n  .  3n 
(**bt)  ^2  *  c  5?  *  d  5? 


where 

a  =■  B  ^1  -  2cte  -  2ct  ^  , 
b  -  2Bcte  , 

c  »  2(cf  +  c  /ir)e8  ,  and 

d  =»  (2/ir)Bc  e  . 

n 


(3) 


The  fundamental  properties  of  the  solutions  of  the  PDE  of  equation  (3) 
depend  on  the  type  of  PDE  (i.e.,  hyperbolic,  parabolic,  elliptic).*  The 
type,  in  turn,  depends  on  the  roots  of  the  radical  as  shown: 

[s2  -  (a+b?)]  1/2  . 

Since  8  is  small,  approximately  1/2,  the  behavior  of  the  solutions  of 
equation  (3)  is  critically  dependent  on  the  sign  of  the  coefficient  of  the 
second-order  space  derivative  expressed  by: 

a+b£  -  8  [l  -  2c^  -  2cte(l-£)]  . 

As  discussed  later,  a+b£  is  typically  less  than  zero  for  most  of  the 
length  of  the  cylinder,  yielding  real  roots  of  the  radical.  Hyperbolic 
systems  of  PDE's  have  real  roots  and  correspond  to  solutions  for  which 
disturbances  on  the  boundaries  and  at  initial  conditions  propagate  along  the 
cylinder . 

A  typical  value  of  2ct£  is  102,  meaning  that  with  no  drogue  or 
instrument  drag  the  PDE  will  be  hyperbolic  at  the  upstream  end  and  remain  so 
over  better  than  99  percent  of  its  length.  The  addition  of  drogue  or 
instrument  drag  could  easily  make  the  PDE  hyperbolic  over  its  entire 
length.  The  physical  significance  of  this  is  that  the  boundary  disturbances 
at  K  “  0  will  propagate  aft  over  most  of  the  cable's  length. 


*A  discussion  of  PDE  types  and  their  characteristics  may  be  found  in 
any  text  that  covers  partial  differential  equations. 
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In  solving  the  equation,  we  make  the  fundamental  assumption  that  the 
solution  of  n(^,T)  is  separable  in  space  and  time,  and  the  time  dependence 
is  harmonic;  that  is, 


n(S,T)  -  v(Oeia)T  , 


(4) 


where  'jo  is  noted  to  be  nond  imensional  and,  in  general,  complex*  and  i 
denotes  the  imaginary  number. 

Substituting  equation  (4)  into  equation  (3),  along  with  the 
transformation  of  the  independent  spatial  variable  (z^  ■  a+b£),  yields  the 
ordinary  differential  equation  (ODE) 


f£2_  dv 
z  dz 


2 

+  b  v 
o 


0  , 


(5) 


where 


2 

a  ■  r-  (c  +  2iu8)  -  1  and 
o  b 


o 


4 <i) 
.2 


(id-'*))  . 


Equation  (5)  is  recognized  as  a  modified  form  of  Bessel's  equation. ^ 


PARTICULAR  SOLUTION 

To  arrive  at  the  particular  solution,  two  required  boundary  conditions 
are  imposed  for  this  second-order  equation.  At  the  upstream  end  of  the 
cylinder  (Jj  *0),  we  impose  the  first  condition  on  the  PDE.  For  £  »  0 
where  0Jo  is  constrained  to  be  real, 

n(?,T>  *  cos  u>oT  • 


Substitution  into  equation  (4)  yields  the  following  for  z^  »  a, 
Re 


Fv(z)1 

cos  a)  T  -  Im  v(z) 

sin  u)  T  *  cos 

L 

0  L  j 

0 

or 


*The  nond  imensional  angular  frequency  of  to  will  be  equated  to  the  real 
nond imens iona 1  frequency  of  0)o  when  we  seek  the  particular  solution  to  a 
forced  vibration  problem.  However,  in  the  complementary  solution  the 
infinite  set  of  Eigen  frequencies  is  necessarily  complex  as  will  be  shown. 
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i 

t 


Re 

Im 


v(z ) 
v  (z) 


1 

0 


Note  that  this  condition  enters  a  forcing  function  through  a  boundary 
condition.  It  is  the  cylinder  response  to  this  forcing  that  is  the  basic 
output  of  this  study  since,  as  discussed  in  the  introduction,  long  cylinders 
are  stable  except,  possibly,  near  the  low  tension  end.  We  will  return  to 
this  point  in  our  discussion  of  the  complementary  solution. 


For  a  boundary  condition  on  the  downstream  end  of  the  cable  (£  *  1),  we 
encounter  a  fundamental  problem  that  is  a  consequence  of  dropping  the 
flexural  rigidity  term  in  the  Paidoussis  momentum  equation.  Since  the 
physical  constraint  at  the  free  end  involves  bending  moments,  the 
second-order  system  (which  allows  no  bending  moments)  is  unable  to 
characterize  the  free  end  physical  condition.  Analytically,  this  problem 
results  from  changing  from  a  fourth-order  to  a  second-order  description.  A 
complete  solution  would  require  matching  the  fourth-  and  second-order 
solutions  in  the  resulting  "boundary  layer"  at  the  free  end.  This  "boundary 
layer"  will  only  apply  to  a  very  short  region  near  the  free  end.  We  will 
not  attempt  this  here,  but  will  seek  another  physical  constraint  which  is 
occasionally  used  in  this  type  of  problem.  Thus  we  will  require  that  the 
lateral  cylinder  displacement  be  bounded;  that  is,  for  0  <  £  <  1, 


j  n(C,x)  <  ®  . 


To  reiterate,  we  seek  the  particular  solution  of  the  ODE.  For 
a  <  z2  <  (a+b), 


i-i  .  ia  p.  .  b2„  .  o 

.2  z  dz  o 
dz 


under  the  conditions  that  for  z2  ■  a, 

v(z)  -  l 

2 

and,  for  a  <  z  <  (a+b),  v(z)  remains  bounded. 


(6a) 


(6b) 


The  range  of  z  given  above  covers  the  entire  length  of  the  cylinder;  that 
is ,  0  <  £  <_  1. 

The  entire  solution  of  v(z)  can  be  written  as  the  sum  of  the  particular 
(v^)  and  complementary  (V2)  solutions;  that  is, 

v(z)  ■  vj(z)  +  V2(z)  .  (7) 

The  complementary  solution  is  addressed  in  the  next  section.  The  development  of 
the  particular  solution  follows. 
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OrCloff  and  Ives  recognized  Chat  a  complete  solution  set  is* 


v(z)  *  zv  Taj  (b  z)  +  BJ  (b  z)l 
L  v  0  -v  O  J 


(8a) 


for  all  values  of  the  complex  order 


1-a 


'  k  (C  +  i28w) 


or 


c  +  (l/ir)c 

_  ,  t  n  .  oj 

v  *  1 - - 


-l 


cte  ^t  cte 


For  all  cases  of  interest  to  this  report,  the  order  is  approximated  as 
real  as  expressed  by: 


(8b) 


which  limits  the  analysis  to  frequencies  much  less  than  ecq/it. 

Bessel  functions  of  complex  order  are  numerically  difficult  to 
evaluate,  and  thus  their  use  in  this  report  would  significantly  increase 
computational  efforts.  Thus  the  region  of  inquiry  is  limited  to  that  for 
which  the  approximation 


c 

V  *  -p  and  p  *  -  >  0  (8c) 


is  appropriate.  (Note  that  the  constraint  required  limits  the  maximum  value 
that  nondimens ional  oj  can  obtain.  This  is  consistent  with  our  previously 
implied  low  frequency  constraint  which  resulted  from  neglecting  the  flexural 
rigidity  term.)  The  solution  may  then  be  written  as 

v(z)  -  z"P  ^AJ 

We  next  address  the  application  of  the  boundedness  constraint  and  the 
upstream  (£  ■  0)  boundary  condition  in  regard  to  assigning  values  to  the  A 
and  B  coefficients  of  equation  (9).  Note  that  z"P  J_p(z)  and 
z~PYp(z)  are  both  singular  td\en  z  is  zero.  Thus,  if  z  is  zero  in  the 


._(b  *) 

-p  o 


BJ  (b  z)l . 
P  o  J 


(9) 


*For  integer  values  of  v,  Jv  must  be  replaced  with  Y^  where  Jv 
and  Y^  are  Bessel  functions  of  the  first  and  second  kind,  respectively. 


region  a  <  z^  <  a  +  b,  Chen  Che  coefficienC  A  must  be  zero  Co  maintain  a 
finice  cransverse  d isplacemenC .  The  region  of  z  in  quescion  is  equivalenC 
Co  Che  region  0  <  £  <  1.  From  equacion  (3)  we  may  wrice 

1  -  2c! 


for  z  *  0. 

Thus  z  cannoC  be  zero  for  0  <  £  <  1  when  2c^  >  1  or  when 
(1  2ct)/2ct£  >  1,  which  implies  chaF  Che  value  of  A  depends  on  Che 

magnitude  of  Che  drogue  or  Cowed  body  Cension.  This  Curns  ouC  noC  Co  be  so 
(i.e.,  A  muse  always  be  zero  for  boundedness),  and  we  may  avoid  Chis 
computational  difficulty  by  reasoning  as  follows:  in  Chis  problem  we  model 
tow-body  effects  as  only  adding  a  constant  drag  to  the  cable  length  of 
interest.  Analytically  (for  our  model)  there  is  no  change  in  either  the  ODE 
or  the  boundary  condition  between  a  cable  length  (L)  plus  a  tow  body  and  a 
cable  length  (L)  plus  an  additional  length  of  the  original  cable  that 
generates  the  same  tension  as  the  tow  body.  These  two  cases  are 
arithmetically  indistinguishable  for  our  model.  It  is,  of  course, 
understood  that  physically  a  tow  body  is  not  completely  described  by  a 
steady  drag  and  that  the  behavior  of  the  two  cases  could  be  different.* 
However,  in  our  current  model  there  is  no  difference.  We  will  also  assert 
that  £ct  »  1  for  practical  cables.  Note  finally,  that  this  discussion 
sheds  light  on  an  issue  left  unresolved  by  Ortloff  and  Ives. 3 

We  will  thus  substitute  a  length  of  the  original  cable  for  the  towed 
package  such  that  the  new  effective  length  (L*)  becomes 


L*  -  L  ♦  d  —  . 

c  ct 


Equation  (3)  then  becomes 


-  8  ^1  -  2ct£  -  2c^  -  6  ^1  -  2cce*^ 


and  equation  (2)  becomes 


E*  -  rx  and  e*  -  L*  . 
L*  3- 


We  now  seek  a  solution  bounded  in  the  region  0  £  <  1.  From  equation 

(3)  we  have 


*This  results  from  the  complex  lift  and  drag  forces  on  the  towed  body, 
which  we  have  not  modeled. 
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1 


2c. 


€* 


for  z  3  0  requiring  that,  because  z  must  be  zero  somewhere  in  the  region  of 
0  <  5*  <  A  must  therefore  be  zero  for  finite  transverse  motion  of  the 
lengthened  cable. 

The  preceding  discussion  justified  making  A  zero  in  general  and  in  no 
way  modified  the  value  of  the  solution  for  0  <  x  <  L.  We  thus  return  to  the 
original  parameter  values  of  equations  (2)  and*  (37* and  note  that  equation 
(9)  becomes 


v(z)  ■  Bz  PJ  (b  z)  .  (10) 

P  o 

Before  evaluating  the  coefficient  B  (a  trivial  matter),  let  us  consider 
the  properties  of  the  above  solution. 

The  values  of  the  solution,  like  the  values  of  a  Bessel  function  with  a 
complex  argument,  are  quite  sensitive  to  the  values  of  the  argument.  We  are 
concerned  with  the  range  and  values  of  the  complex  quantity,  b0z .  From 
equation  (5b)  we  may  write 


b  z  3  + 
o 


*4/2 

-  TT72 


J.  d_^o 
V  1  b  b 


(11a) 


Substituting  parameter  values  for  the  b  and  d  coefficients  from 
equation  (3b)  yields 


2  ml 


1/2 


V  -  -  ?7T 


VC  (1) 

.  n  o 

1  TTCt  '  20ct, 


To  the  same  level  of  approximation  as  was  required  to  make  the  Bessel 
function  index  real,  we  have  from  equation  (8b) 


b  z  s 

o 


4uopi 


which  becomes 


b  z 
o 


i(3ir/4) 

e  z 


for  z2  <  o 


(lib) 
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.#* 

if 

$ 

•V.r 


•.ft 


:y; 

y/ 

vr 


for  >  o.  Note  that  the  choice  of  sign  in  equation  (11)  is  unimportant 
because  of  the  synmetry  of  the  basic  ODE  in  equation  (5)  to  b0  and  z. 

Figure  2  is  a  plot  of  the  real  and  imaginary  values  of  the  argument, 
bQz,  in  equation  (11a)  as  the  frequency  (w0)  and  space  (£)  variables  are 
changed.  The  solid  lines  illustrate  the  values  of  b0z  at  the  upstream 
(£  *  0)  and  downstream  (K  ■  1)  boundaries  as  (,)0  is  increased.  The  dotted 
line  is  the  locus  of  values  of  b0z  for  a  constant  u0  and  C  changing  from 
0  to  1.  Note  that  the  figure  is  for  the  case  in  which 


for  S  *  0,  and 


z2  >  0 


for  K  ■  1.  The  range  of  values  for  bQz,  which  is  consistent  with  the 
assertion  of  equation  (8b)  that  is,  equation  (lib),  is  shown  by  the 
straight  line  portion  of  the  solid  line.  Also  shown  in  figure  2  is  the 
position  of  the  system  poles  (Eigenvalues)  as  derived  in  the  next  section  of 
this  report.  For  the  range  of  wo  for  which  equation  (lib)  applies,  the 
argument  bQz  is  well  removed  from  any  resonant  effects  as  is  seen  by  the 
separation  of  the  values  of  b0z  and  the  system  poles  in  figure  2. 

Evaluation  of  equation  (10)  for  the  frequency  range  of  w0  for  which 
equation  (lib)  applies,  results  in  shapes  v(z)  that  always  oscillate  and 
decay  as  5  goes  from  0  to  -a/b,  if  -a/b  is  less  than  one,  and  from  0  to  1, 
if  -a/b  is  greater  than  one.  For  typical  parameter  values  of  interest,  -a/b 
is  approximately  equal  to  or  greater  than  unity. 

It  is  instructive  to  analyze  the  small  argument  case  (bDlzl  <<:  1), 
which  would  occur  at  small  frequencies  of  u  <K  d,  by  evaluating  equation 
(10)  using  the  first  term  of  the  defining  power  series  for  a  Bessel 
function. ™  This  gives  the  equation 


i(p1T/4) 


bpT(p+l) 


for  0  £5^1,  where  ^  (  )  denotes  a  gamma  function.  Equation  (12)  simply 

states  that  at  a  very  low  frequency  the  entire  cylinder  duplicates,  without 
delay,  the  motion  of  the  tow  point  (5  ■  0).  Incorporating  the  second  term 
of  the  defining  power  series  yields 


i(V4) 


bpF(  p+i) 


[*  '  1  ^  <**b5)] 


The  term  in  brackets  may  be  approximated  as 

exp  T-i(*.  +  $„5)1  , 
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where 


<t> 


1 


2 


o>ad 

b2(p+l) 

cod 

b(  p+1 )  * 


Substituting  the  values  of  equations  (3b)  and  (8c)  into  the  above  expression 
for  <$2  shows  that  equation  (13)  represents  a  one-dimensional  wave  propa¬ 
gating  downstream  with  no  attenuation  and  a  nondimensional  phase  speed  of 


(c  /ttc.  )  +  1 
n _ c _ 


(14) 


This  is  the  "water  pulley"  effect  described  by  Kennedy. ^ 

Now  let  us  complete  the  particular  solution  of  equation  set  (6)  by 
evaluating  the  complex  constant  B  from  equation  (10)  using  the  boundary 
condition  of  equation  (6b). 


B 


aP/2 

J  (b  a172) 
P  o 


(15) 


which  allows  the  particular  solution  to  be  written  as  follows.  For 
a  <  z?-  <  a+b, 


Vj  (z) 


J0>  z) 

P 

T  l/2 
J  (b  a 
p  o 


(16) 


COMPLEMENTARY  SOLUTION 

The  complementary  solution  results  from  solving  equation  (6a)  with  the 
boundary  conditions 

v2(z)  »  0  (17a) 

for  z^  •  a,  and  V2(z)  remains  bounded  for 

a  <  z2  a+b  .  (17b) 

Equation  (10)  is  the  solution.  Thus,  with  the  upstream  boundary  condition 
in  equation  (17),  we  require  that  b0a*/2  be  equal  to  the  zeros  of  the 
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$ 

4 

w* 


i 


& 


>•1 


V* 


Bessel  function  of  order  p.  Since  Bessel  functions  with  a  positive  real  index 
are  real, 12  we  may  write  this  equation  as  follows.  For  n  ■  1,  2,  .  .  .  , 


where  5n  are  the  infinite  real  zeros  of  JpCboal/2).  The  complex 
frequencies  required  to  meet  the  boundary  condition  of  equation  (17)  may  be 
found  by  substituting  equation  (5)  into  equation  (18)  and  solving  for  the 
Eigen  frequencies  (u^)  as  follows: 


id  .  1  (  .1  b2  .2\ 

r  -i  p  -rV 


Let  us  examine  the  Eigen  frequencies  by  substituting  equation  (3b)  into 
equation  (19).  After  some  manipulation,  the  Eigen  frequencies  may  be 
written  as 


“n  ‘  1  TT2  [‘I  (> 


<V»> 


2  .  1/2 


)  ] ' 


(20a) 


where  D  is  a  drag  term  defined  as 


D  -  2c„e  +  2c„  . 
t  c 


Only  the  Eigen  frequencies  that  meet  the  requirement  of  equation  (8b) 
can  be  used ,  that  is , 


Thus,  for  each  U)_  we  must  test  the  following  condition 


i8p  £  1  ♦  ^ 


(yP> 


2  v  1/2 


S(l-D) 


(20b) 


This  is  a  very  restrictive  condition.  To  illustrate  this  point  we  will 
examine  the  condition  for  the  first  (smallest)  zero  (6j).  Examination  of 
equation  (20b)  shows  that  if  does  not  meet  the  condition  no  other  zero 
will  either. 


The  first  zero  lies  in  the  interval!^ 


p  <  <  ^2(p+l )(p+3)J 
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For  large  p,  5^  *  >Jl  p.  So,  within  a  factor  of  two  we  may  write 
6^  -  p,  and  equation  (20b)  becomes 


I  is"  [l-  (l  *  3TT5t)1/2  ]  I 


«  p 


A  systematic  evaluation  of  the  above  equation  reveals  that  only  when 
S(l-D)  »  1  does  the  Eigen  frequency  meet  the  equation  (8b)  condition.* 
This  statement  may  be  easily  generalized  to  read: 


0) 


n 


ic8ct/* 


U„/P>2 

B(D-l) 


for 


B(D-l)  " 


only. 


This  may  be  simplified  to  read: 


.  £Cn  (5n/p)2 

Wn  "  1  2tt  Td^D 


(21) 


(«n/P)2 

for  a To^TT  k<  1  • 


To  sunmarize,  all  of  the  Eigen  frequencies,  which  we  can  calculate, 
indicate  that  the  system  will  be  temporally  stable  since  the  solution 
consists  of  a  sum  of  terms  in  the  form 

B  z'P  J  (5  )e%T  ,  (22) 

n  p  n 


where  a>Q  is  a  positive  imaginary  number. 


Of  interest  to  our  present  study  is  the  damping  time  associated  with 
each  Eigen  frequency  because  this  would  be  the  time  required  for  ship 
maneuvering  effects  to  have  died  out  after  assuming  a  constant  average 
heading.  It  is  useful  to  define  a  damping  time  (t0)  as 


Im 


1  . 


*Note  that  the  fact  that  3  is  approximately  1/2  figures  into  this 
conclusion. 
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The  slowest  mode  (i.e.,  longest  time  constant)  for  T  >1,  is  given  by 


1 


Tmax  min  [ Im(Wn)]  ’ 


(23) 


where  min  [imCa^)]  »  Im(oo].)* 
approximated  as  before  as  6^  : 


The  first  zero  (5^)  of  Jp  may  be 
p.  From  equation  (22)  we  have 


In  real  (dimensional)  time 


^  To \nax 


max 


2tt(D-1) 

ec 


(24) 


2„<D-1>L  2’'(°-I)d= 


Uc 


(25) 


Typical  values  for  these  parameters  are  ct  ■  .0025,  cn 
e  *  75000,  dc  *  .01m,  U  ■  3m/sec;  therefore, 

(t  )  2:  150  seconds  . 

o  max 


.07,  ct  -  57, 


Ortloff  and  Ives^  concluded,  wrongly,  that  the  cylinder  would  always 
be  temporally  unstable.  Lee  and  D'Appolito^  corrected  Ortloff' s  and 
Ives'  error  and  concluded  that  if  some  minimum  tension  were  exceeded  then 
the  system  would  be  temporally  stable.  However,  neither  of  these  papers 
tested  the  Eigen  frequencies  for  conformance  with  the  condition  of  equation 
(8b).  The  above  discussion  shows  that  the  region  identified  by  Lee  and 
D'Appolito  (D  <  1)  as  being  unstable  is  not  necessarily  so  because  for  that 
range  of  tension  the  order  of  the  Bessel  function  cannot  be  approximated  as 
real.  Furthermore,  some  of  the  Eigen  frequencies  are  not  computable  because 
the  Bessel  function  index  is  required  to  be  real,  so  no  definite  temporal 
stability  statement  is  possible.  The  point  appears  to  be  academic,  however, 
because  there  is  no  experimental  reason  to  believe  that  an  instability 
exists. 


ANALYTICAL  EXAMPLES 

A  qualitative  display  of  cylinder  dynamics  is  shown  in  figures  3 
through  5.  These  figures  were  obtained  by  numerically  evaluating  equations 
(4)  and  (16)  for  a  range  of  values  for  £  and  j .  The  displays  show  time 
histories  of  the  particular  solution  (i.e.,  the  cylinder  transverse  shape  at 
three  different  frequencies).  In  all  cases,  the  sinusoidal  boundary 
condition  imposed  on  the  tow  point  propagates  attenuatively  back  through  the 
cylinder.  As  expected,  the  attenuation  rate  clearly  increases  with  the 
frequency. 
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Let  us  define  the  transfer  function  [h(oj)]  to  -be  the  complex  ratio  of 
the  Fourier  Transforms  of  the  crosstrack  displacements  at  an  arbitrary  point 
on  the  cylinder  (£)  to  the  upstream  boundary  (£  *  0).  Using  equation  (4) 
this  ratio  mav  be  written  as 


H(u>)  = 


where,  from  equation  (6b),  vi(0,ujo)  =  1.  Thus, 


H(u))  -  |Vl(5,U)o)|  e^’V  , 


where 


•i  fI"(vi>l 


The  transfer  function  (ordinate  of  figure  6)  is  the  ratio  of  the 
amplitude  of  the  sinusoidal  motion  of  the  cylinder  observed  at  some  spatial 
point  (£)  to  that  imposed  on  the  cylinder  at  the  forward  boundary  (£  *  0) 
for  a  given  frequency  (u)0).  From  figure  6  we  note  that  the  observer  at  £ 
sees  the  transfer  function  diminish  (from  unity)  as  frequency  increases 
(i.e.,  the  cylinder  passes  low  frequencies  and  rejects  high  frequencies). 

A  point  of  constant  phase  obeys  U)T  -  \|>(£)  *  constant.  Thus  that  point 
will  have  a  velocity  given  by  0JdT  -  t^'dC  *  0,  and  the  phase  velocity  (c)  is 


d£  a  u) 
dx  r  * 


(27a) 


Unfortunately,  the  spatial  derivative  (^’)  could  not  be  measured 
directly  in  the  companion  experiment  described  in  the  next  section.  Thus  we 
define  c*  as  the  phase  velocity  measured  between  £  and  £  equal  to  zero,  as 
expressed  by: 

c*  -  SS.  .  (27b) 


The  c*  is  a  function  of  £  and  u,  requiring  a  three-dimensional  display 
of  the  nondimens ional  phase  velocity  as  given  in  figure  7.  Because  of  the 
difficulty  of  displaying  this  function,  two  different  viewing  angles  of  the 
same  function  are  shown  in  the  figure.  At  low  frequencies  the 
nond imens iona 1  phase  speed  is  given  by  equation  (14)  as  c*  ■  1  +  p“l. 
Typical  values  of  p  are  greater  than  10,  which  makes  c*  nearly  one  and 
independent  of  5  as  was  predicted  by  equation  (13).  As  the  frequency 
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Figure  6.  Space  Frequency  Plot  of  Transfer  Function 


increases,  c*  increases  (dispersive ly)  with  larger  increases  occurring 
nearer  5  «  0  than  £  »  1.  This  picture  is  totally  consistent  with  the 
physical  description  given  by  Kennedy^  and  sunmarized  as  follows:  at  low 
frequencies  the  hydrodynamic  forces  dominate  over  the  cylinder- fluid 
inertial  and  cylinder  tension  stiffness  terms,  causing  the  disturbances  to 
propagate  downstream  at  the  flaw  velocity  (i.e.,  c  •  1).  Because  tension 
stiffness  is  negligible,  the  phase  velocity  is  independent  of  tension  and 
thus  independent  of  5.  As  frequency  increases,  both  the  above  mentioned 
forces  (inertial  and  tension  stiffness)  become  important  and  the  cylinder 
responds  with  "string-like"  waves  having  propagation  velocities  greater  than 
flow  velocity.  The  velocity  of  these  waves  is  proportional  to  the  square 
root  of  the  cylinder  tension  and  thus  decreases  as  £  goes  from  0  to  1. 


HORIZONTAL-VERTICAL  INTERACTION 

Up  to  this  point  in  the  report  we  have  considered  the  cable  of  interest 
to  be  a  neutrally  buoyant  flexible  cylinder.  Actually,  most  cables  are 
negatively  buoyant.  Rather  than  undertake  the  monumental  task  of  rederiving 
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and  then  solving  equation  (1)  in  three  dimensions,  we  propose  to  modify  the 
coefficients  of  equation  (l)  to  account  for  the  new  forces  associated  with 
the  vertical  angle  assumed  by  a  negatively  buoyant  cable.* 

Figure  8  illustrates  a  segment  of  a  cable  being  towed  in  the  negative  y 
direction,  having  a  horizontal  angle  with  the  flow  of  9,  and  having  a 
depression  angle  <J>  due  to  uniform  cable  weight.  The  base  vector  in  the 
cable  direction  is  c  and  may  be  shown  as: 

-  A  A  A  .  , 

c  *  i  sin  9  cos  $  +  j  cos  9  cos  <j>  +  k  sm  <j>  ,  C 28 J 

AAA 

where  i,  j,  and  k  are  base  vectors  in  a  Cartesian  coordinate  system.  The 
fluid  flow  vector  is 

U  *  'jU  .  (29) 


FLOW  DIRECTION 


BASE  VECTOR  DIAGRAM 


Figure  8.  Horizontal-Vertical  Force  Diagram 

*This  approach  was  originally  brought  to  the  author's  attention  by 
Dr.  R.  Kronauer  of  Harvard  University. 
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The  tangential  and  normal  components  of  velocity  relative  to  the  cable  are 
found  from  the  following  defining  equations: 


UJ.  -  (T  *  £) c 


U 


N 


(30) 


After  some  tedious  but  straightforward  algebra  we  may  express  the  normal 
component  of  velocity  as 


-  ^  cos  9  cos^  $  sin  9 

+  ^(1  -  cos^  9  cos2  <J>) 

-  1c  cos  9  cos  ♦  sin  $ 


The  magnitude  squared  value  is 

„  -2 

N 


2  2 

■  1  -  cos  9  COs  $  , 


(31) 


(32) 


Equations  (31)  and  (32)  show  that  the  component  of  flow  normal  to  the 
cable  is  a  function  of  both  the  horizontal  and  vertical  angles.  Thus,  the  x 
component  of  the  normal  flow,  which  is  what  interests  us,  is  also  a  function 
of  both  angles.  This  means  that  there  is  a  horizontal  force  on  the  cable  as 
a  consequence  of  the  vertical  angle  4>.  We  now  proceed  to  calculate  that 
force.  The  normal  force  vector  is 

"fJ,  -  R(1  -  cos2  9  cos2  <)>)n  ,  (33) 

where  n  is  the  base  vector  normal  to  the  cable  (i.e.,  n  *  ^c  *  0)  and 

R  -  I  Pc  d  U2  .  (34) 

2  n  c 


The  quantity  n  is  given  by 


A 

n 


0(1  -  cos2  9  cos2  $)*^2 


(35) 
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To  find  the  horizontal  component  of  force,  substitute  equation  (35) 
into  equation  (33)  and  identify  the  x  component  of  the  result.  Again,  after 
some  algebra,  we  have 


2  22  1/2 
(F„)  *  -R  cos  0  cos  *  sin  9  (1  -  cos  9  cos  4>)  , 

N  x 


(36) 


which  compares  with  the  neutrally  buoyant  case  ($  *  0)  of 

(Vx  I  $-0  *  "R  C09  9  si°2  9  *  (37) 

(Ffj)x  J  c^e  normal  hydrodynamic  drag  term  of  equation  (1). 

A  comparison  of  equations  (36)  and  (37)  shows  that  increasing  the 
depression  angle  <J>  to  about  40  degrees  increases  the  horizontal  force  by  the 
amplification  factor 


2  ...  2  a  2  , .  1/2 
cos  4>(1  -  cos  9  cos  q>) 

sin  9  ’ 


(38) 


where 


<F»>, '  • 

A  convenient  way  to  incorporate  the  amplification  of  equation  (38)  into 
the  solution  of  equation  (5)  is  to  increase  the  coefficient  of  normal  force; 
that  is,  c£  will  be  used  in  place  of  cn  where 


* 

c  -  Ac 
n  n 


(39) 


Note  that  singular  behavior  of  A  for  zero  6  has  no  physical 
significance  since  (Fjj)x  is  zero  at  zero  9.  This  is  addressed  in 
appendix  B. 

Physically  we  have  shown  in  this  section  that  the  normal  hydrodynamic 
force  on  the  cable,  resulting  from  a  vertical  depression  angle,  will  have  a 
horizontal  component  when  the  cable  is  angled  to  the  flow  in  the  horizontal. 
This  describes  a  horizontal-vertical  interaction  of  forces.  The  vertical 
cable  angle  establishes  a  flow  across  the  cable  that  results  in  a  force 
aligned  with  the  flow  in  the  absence  of  a  horizontal  inclination  of  the 
cable.  However,  if  the  crosstrack  meandering  of  the  towing  vessel  causes  a 
horizontal  attack  angle  on  the  cable,  then  the  normal  force,  previously 
aligned  with  the  flow,  is  "tipped"  in  the  horizontal  direction.  Equation 
(39)  modifies  the  normal  drag  coefficient  to  account  for  this  additional 
horizontal  force. 
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EXPERIMENTAL  PROGRAM 


DESCRIPTION 


During  Che  experimental  program,  a  thin,  16-am  diameter,  steel  cable 
was  towed  from  a  surface  vessel,  and  the  vessel  was  maneuvered  to  force 
sinusoidal  horizontal  transverse  motion  at  the  upstream  end  of  the  cable. 
During  each  run,  the  cable  was  towed  through  a  number  of  cycles  (at  least 
four).  Each  run  was  defined  by  a  set  of  three  parameters  varied  in  a 
factorial  type  experiment  design. ^  The  variables  were  ship  speed,  the 
length  of  cable  between  the  ship  and  the  pinger,  and  the  frequency  of  the 
oscillation.  These  experimental  parameters  are  illustrated  in  figure  9. 


The  transverse  motions  of  both  the  ship  and  a  point  on  the  cylinder 
were  carefully  measured  by  acoustically  tracking  pingers  attached  to  the 
ship  and  the  cable  at  the  point  to  be  measured.  As  shown  in  figure  10,  the 
tracking  was  done  from  spatially  fixed  hydrophones.  The  ship  maneuvers 
simulated  a  sinusoidal  upstream  boundary  condition,  and  the  response  of  the 
system  (to  the  stimulation)  was  measured  by  the  track  of  the  pinger  at  the 
end  of  the  cable.  A  neutrally  buoyant  cylinder  supplied  tension  (drogue)  to 
the  end  of  the  cable. 


The  objective  of  the  data  analysis  was  to  compare  the  experimental 
results  and  the  theoretical  results  discussed  in  this  report.  An  exact 
comparison  was  possible  for  the  magnitude  of  the  transfer  function,  equation 
(26),  and  the  phase  velocity,  equation  (27).  Experimentally,  this  was 
determined  by  calculating  the  Fourier  transform  of  the  time  series  of  the 
crosstrack  displacements  (see  figure  9)  for  both  pingers.  The  transfer 
function  magnitude  |h(w)  |  and  phase  speed  (c)  were  then  determined  as  in 
equations  (26)  and  (27b)  for  the  single  frequency  at  which  the  system  was 
forced . 


As  shown  in  figure  10,  the  thin,  flexible  cable  was  negatively  buoyant 
as  a  consequence  of  the  uniform  linear  weight  of  the  steel  cable.  Because 
of  the  uniform  linear  weight,  the  cable  tows  with  a  constant  depression 
angle  4*  (called  the  "critical  angle")  from  the  surface  down  to  a  depth  at 
wh ich  the  drogue  tension  causes  the  depression  angle  to  transition  from  $  to 
zero.  Knowing  the  cable  length*  and  the  pinger  depth  from  a  depth  sensor 
(see  figure  10)  allows  an  approximate  calculation  of  the  average  depression 
angle  as  follows 


*  =>  5in-l  f~ pinger  depth'  _ 

L  c/s  J 


This  depression  angle  varied  from  4  to  9  degrees  during  the  experiment.^ 


*It  is  important  to  differentiate  between  the  cable  length  (Lc/s)  of 
figure  10  and  the  horizontal  projection  of  cable  length  (L)  of  figure  9.  It 
is  the  length  of  the  cable  in  the  direction  of  fluid  flow  (L)  which  appears 
in  equation  (1). 
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When  the  experiment  was  performed,  the  Rispin^  results  were  not 
available;  consequently,  the  horizontal  angle  limitations  on  the  Paidoussis 
linearization  were  not  known  (refer  to  discussion  in  introduction'. 
Therefore,  the  angle  0  frequently  exceeds  the  2.5  to  3  degree  threshold 
found  by  Rispin.  Note  from  figure  9  that  0  will  vary  sinusoidally  with  a 
maximum  value  of  ka  and  an  average  magnitude  of  2/-rrka.  Experimentally  the 
maximum  angles  varied  from  1.8  to  29.8  degrees,  and  the  average  magnitudes 
varied  from  1.1  to  18.9  degrees.  From  equation  (38),  it  is  seen  that  the 
amplification  of  the  horizontal  force  applied  to  the  cable  is  a  function  of 
both  0  and  <j>.  Figure  11  illustrates  the  range  of  amplification  that  occurs 
for  the  experimental  range  of  0  and  ^  cited  above.  Because  0  varies  in  both 
space  and  time,  so  will  the  amplification  factor.  An  effective  value  will 
be  sought  experimentally  rather  than  analytically. 

Quantitative  statements  concerning  the  experimental  accuracy  will  not 
be  made  for  the  following  two  reasons.  The  tracking  accuracies  of  the 
acoustic  tracking  range  are  not  well  understood;  this  is  particularly  true 
regarding  how  those  tracking  errors  propagate  through  the  data  analysis 
used.  Tracking  errors  represent  noise  in  the  process  that  is  not  known.  In 
general,  estimates  of  a  transfer  function  based  on  a  finite  length  time 
series  will  contain  errors  even  when  the  time  series  contains  no  noise. ^ 
This  error  is  a  function  of  both  the  transfer  function  to  be  measured  and 
the  test  signal  used  in  the  measurement.  Quantitative  evaluation  of  this 
error  is  not  trivial  and  was  not  undertaken.  However,  it  is  thought  that 
neither  of  these  errors  will  bias  the  transfer  function  estimates  and  that 
on  the  average  these  errors  should  contribute  only  random  (zero  mean)  error. 

Before  an  analytical  comparison  can  be  made  with  the  data,  the  "drogue" 
factor  (Ct)  needs  to  be  evaluated  from  the  following  expression 

c r 

t  ^  drag  of  drogue 
cte  drag  of  cable 

The  Ct>  from  equation  (1),  accounts  for  the  component  of  the  total 
drag  that  is  constant  over  cable  length  L;  i.e.,  ct  is  proportional  to 
the  component  of  tension  that  varies  linearly  over  L,  and  ct  is 
proportional  to  the  component  that  is  constant  over  L.  The  drags  of  the 
drogue  and  cable  were  calculated  using  the  appropriate  diameters,  lengths, 
and  standard  formulas,  i.e.,  equation  (lb). 


RESULTS 

The  results  of  the  data  analysis  are  shown  in  figures  12  and  13. 

Tables  A-l  and  A-2  of  appendix  A  summarize  the  condition  of  each  one  of  the 
experimental  points  for  each  of  the  two  experiments.  The  large  dots  on 
figure  12  are  the  experimental  values  of  the  transfer  function  magnitude, 
previously  described,  plotted  as  a  function  of  nondimens ional  frequency. 

The  dots  in  figure  13  show  the  phase  velocity,  nond imens ional ized  by  the 
towing  speed  through  the  water,  also  plotted  as  a  function  of  the 
nondimensional  frequency. 
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The  solid  lines  on  figures  12  and  13  are  theoretical  curves  resulting 
from  equations  (16)  and  (27b).  All  of  the  parameters  used  in  the 
calculation  are  reasonably  well  established  from  previous  experience  except 
for  the  normal  drag  coefficient  (cn),  which  is  the  coefficient  involved  in 
the  linearized  normal  drag  term.  It  is  the  value  of  this  term  that  is  least 
certain.  The  Rispin^  experiment  showed  that  for  shallow  attack  angles, 
less  than  2.5  to  3  degrees,  the  form  assumed  by  Paidoussis  was  correct  and 
had  a  coefficient  value  of  cn  equal  to  .070.*  From  equation  (39)  and 
figure  11  we  would  anticipate  that  the  effective  coefficient  (c^) 

(which  accounts  for  the  effect  of  a  steady  vertical  tow  angle),  would  be  two 
to  three  times  the  Rispin  coefficient.  To  test  this,  we  next  determined  the 
value  of  c„  that  yielded  the  best  visual  fit  to  both  the  amplitude  and 
phase  data  (figures  12  and  13)  simultaneously.  The  value  found  was  0.16 
compared  with  Rispin’ s  value  (.070).  From  equation  (39),  c^  a  0.16 
implies  an  effective  amplification  value  of  2.29.  This  value  falls  within 
the  independently  predicted  value  range  shown  in  figure  11.  The  theoretical 
curves  obtained  this  way  are  labeled  "Rispin  Normal  Drag  Coefficient 
Modified  for  Cylinder  Depression  Angle"  in  figures  12  and  13.  The 
sensitivity  of  the  complex  transfer  function  to  changes  in  the  value  of 
c*  was  examined  during  the  above  curve  fitting  procedure.  It  was  found 
that  factor  of  two  changes  in  c^  yielded  changes  in  { H ( to)  j  and  c*/U  of 
less  than  +0.1  in  both  cases. 

The  rest  of  the  analytical  coefficients  are  £  *  75000,  cn  *  .2020, 

£  ■  1.0,  8  ■  0.5,  ct  ■  .0101,  and  cj.  =  228. 

Other  theoretical  curves  were  calculated  using  two  different  values  of 
c^i  one  curve  represents  the  best  visual  fit  to  the  transfer  function 
magnitude,  and  the  other  represents  the  best  visual  fit  to  the 
nondimens ional  phase  velocity.  These  are  labeled  "Normal  Drag  Coefficient 
for  Best  Data  Fit"  in  figures  12  and  13. 

Figure  12  shows  that  use  of  the  Rispin  parameter  overestimates  the 
measured  transfer  function  magnitude  and  that  a  smaller  value  of  cn  is 
necessary  to  match  the  data.  Since  the  only  attenuative  mechanism  available 
is  the  normal  hydrodynamic  drag,  which  introduces  a  velocity  damping,  this 
result  would  appear  to  run  counter  to  intuition.  That  is,  attenuation 
increases  with  an  increase  in  normal  drag.  To  see  the  source  of  this 
apparent  contradiction,  let  us  recall  that  the  normal  drag  force  is  modeled 
as 

(*♦«£) 

in  equation  (1),  where  the  flow  across  the  cable  is  3y/ 3t  +  U(3y/9x). 
Presumably  an  increase  in  cn  should  increase  the  attenuation  for  a  given 
frequency.  Let  us  examine  the  normal  force  when  y  is  a  wave  traveling 
downstream  with  the  phase  velocity  c(^),  as  expressed  by: 


*This  value  will  be  referred  to  as  the  Rispin  coefficient. 
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where  k  *  °>  -r  . 

c(oj) 


y  *  Yei(a,C"kx> 


Then, 


F»  -  v1  f  (f-1)  ei<“t"tot)  • 

Thus  we  see  Chat  Che  phase  velocity  must  be  different  than  U  for  a  damping 
force  to  exist.  But  we  recall  from  the  discussion  of  equation  (14)  that  the 
phase  velocity  is  proportional  to  1  +  P~l.  A  decrease  in  en  is  thus 
needed  to  increase  the  attenuation.  Likewise,  in  figure  13,  an  increase  in 
cn  was  necessary  to  best  fit  the  phase  velocity  data.  The  absence  of  a 
single  value  of  cn  that  best  fits  both  the  amplitude  and  the  phase  of  the 
measured  transfer  function  illustrates  the  difficulty  of  the  linear  model  in 
matching  data. 

The  data  was  next  examined  to  determine  whether  a  better  fit  to  the 
curves  calculated  using  the  modified  Rispin  coefficient  would  be  possible  by 
editing  out  those  runs  that  had  very  large  8.  The  results  were  negative. 
Other  attempts  at  observing  a  nonlinear  amplitude  effect  have  yielded 
similar  negative  correlations.  Presumably  a  better  fit  to  theory  would  be 
possible  with  a  more  rigorous  handling  of  the  three-dimensional  dynamics  and 
a  more  exact  treatment  of  the  nonlinear  hydrodynamic  normal  drag.  However, 
the  current  theory  does  agree  with  experimental  values  to  roughly  25  percent. 


SUMMARY 


A  linear  theory  describing  the  transverse  motion  of  a  cable  being  towed 
by  a  vessel  imperfectly  attempting  to  steer  a  straight-line  course  was 
studied  both  theoretically  and  experimentally.  The  analytical  approach  was 
to:  (1)  expand  the  existing  theory  of  the  transverse  motion  of  neutrally 

buoyant,  long,  flexible  cylinders  in  a  uniform  axial  flow,  and  (2)  modify 
the  coefficients  of  this  theory  to  account  for  the  horizontal-vertical 
interaction  of  forces -resulting  from  the  steady  vertical  towing  angle 
associated  with  the  towing  of  negatively  buoyant  cables.  This  approach  was 
taken  to  avoid  the  difficult  problem  of  deriving  and  solving  the  appropriate 
governing  equation  in  three  dimensions. 

The  results  of  the  numerical  calculations  show  that  for  parameters 
typical  of  long  cables  towed  at  sea,  disturbances  occurring  on  the  upstream 
boundary  (tew  point)  propagate  dispersively  and  attenuatively  downstream. 
Both  the  attenuation  and  the  phase  velocity  of  the  disturbance  are  functions 
of  the  frequency  content  of  the  disturbance  and  the  distance  between  the  tow 
point  and  an  observer.  These  examples  served  to  quantify  a  physical 
description  of  the  process  previously  given  by  Kennedy.  H  In  that 
description  it  was  shown  that  at  low  frequencies,  where  inertial  and  tension 
forces  are  dominated  by  hydrodynamic  drag  terms,  the  attenuation  is  small 
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and  the  phase  velocity  is  near  the  fluid  flow  velocity.  Both  quantities 
change  little  spatially  because  the  hydrodynamic  forces  are  spatially 
uniform.  As  the  frequency  increases  enough  for  the  inertial  forces  and 
tension  forces  (because  of  the  increased  curvature  of  the  cylinder)  to 
become  significant,  the  attenuation  and  phase  velocities  increase  and  both 
become  spatially  dependent  as  a  consequence  of  the  linearly  changing 
tension.  The  theory  developed  in  this  report  shows  that  there  are  no 
resonances  encountered  by  excitation  by  real  frequencies. 

The  temporal  stability  of  these  cylinders  was  also  investigated  with 
limited  results.  A  basic  low  frequency  assumption  is  made  throughout  this 
report  which  limits  one's  inspection  of  all  Eigenvalues  of  the  system. 
However,  those  Eigenvalues  that  lead  to  Eigen  frequencies  within  the  region 
for  which  the  present  theory  applies  indicate  a  temporally  stable  system. 
This  is  thought,  in  general,  to  be  true  because  no  instabilities  have  been 
experimentally  observed. 

It  was  shown  in  this  study  that  if  a  cylinder  is  statically  inclined 
vertically  as  well  as  dynamically  inclined  horizontally,  there  is  an 
effective  amplification  of  the  horizontal  component  of  normal  force  on  the 
cable.  This  accentuates  the  hydrodynamic  forces  and  thus  causes  the  "water 
pulley"  effect. 


An  experiment  was  performed  at  sea  to  test  the  theory.  The  parameter 
values  used  in  the  calculations  were  standard  except  for  the  hydrodynamic 
normal  drag  term  coefficient.  Ultimately,  the  value  for  the  coefficient 
measured  in  the  laboratory  by  Rispin^,  amplified  to  account  for  the 
horizontal-vertical  interaction  of  forces  discussed  above,  yielded  a 
solution  that  fit  the  data  to  within  25  percent.  We  are  therefore  confident 
that  we  had  modeled  the  dominant  physical  mechanisms.  Better  individual 
fits  to  the  magnitude  and  phase  components  of  the  transfer  function  were 
possible  by  changing  cn  for  each  component.  The  absence  of  a  single  value 
of  cn  that  best  fit  both  the  amplitude  and  phase  of  the  measured  transfer 
function  illustrates  the  difficulty  of  the  linear  model  in  matching  data.  A 
better  fit  to  theory  would  probably  be  possible  with  a  more  rigorous 
handling  of  the  three-dimensional  dynamics  and  a  more  exact  treatment  of  the 
nonlinear  hydrodynamic  normal  drag. 
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APPENDIX  A 

SUMMARY  OF  EXPERIMENTAL  PARAMETERS  AND  RESULTS 


Tables  A-l  and  A-2  sumaarize  Che  parameters  and  results  of  two 
experiments.  Experiments  1  and  2  were  performed  in  April  1979  and  May  1980, 
respectively.  The  following  symbols  are  used  in  the  tables: 


f 

L 

a 

X 

U 


is  forced  vibration  frequency  in  mHz, 

is  horizontal  distance  between  acoustic  pingers  in  meters, 

is  the  peak  amplitude  of  the  forced  vibration  in  meters, 

is  the  wavelength  of  the  forced  vibration  in  meters, 

is  the  average  velocity  of  the  water  relative  to  the  cable  in 
meters  per  second , 


c*  is  phase  velocity  measured  between  the  acoustic  pingers 
nondimens ionalized  by  flow  velocity,  and 

|H|  is  the  magnitude  of  the  transfer  function  between  the  two 
sensors . 


The  major  objective  of  experiment  1,  the  initial  test,  was  to  determine 
the  feasibility  of  the  experimental  concept  and  to  develop  some  experimental 
procedures.  During  experiment  1  there  was  only  one  tow  vessel  speed, 

U  »  3m/s,  and  no  modifications  of  that  number  were  applied  to  account  for 
vertical  current  shear.  By  the  time  the  data  analysis  for  the  second 
experiment  took  place,  more  accurate  ways  of  determining  U  and  L  were 
ava  ilable. 


Table  A-l.  Experimental  Parameters  and  Results 
for  Experiment  No.  1 


Run 

No. 

f 

(mHz) 

L 

(m) 

H 

c 

a/X 

1 

0.5 

1135 

0.995 

1.042 

0.02 

2 

1.0 

1135 

0.917 

1.034 

0.04 

3 

2.0 

1135 

0.681 

1.046 

0.07 

4 

0.5 

635 

0.830 

0.898 

0.02 

5 

1.0 

635 

0.736 

1.094 

0.04 

6 

2.0 

635 

0.414 

1.083 

0.07 

A-l 
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Table  A-2.  Experimental  Parameters  and  Results 
for  Experiment  No.  2 


Run 

No. 

L 

(m) 

a/X 

U 

(m/ s ) 

f 

(mHz) 

H 

c* 

1 

1186.9 

0.0097 

3.10 

0.97 

0.8223 

1.00 

2 

1186.4 

0.0195 

3.03 

2.07 

0.3804 

1.10 

3 

1186.0 

0.0395 

2.98 

3.79 

0.1075 

1.05 

4 

1185.7 

0.0829 

2.94 

6.84 

0.0085 

1.18 

5 

1193.8 

0.0065 

5.10 

0.37 

0.9671 

1.19 

5A 

1194.1 

0.0065 

5.25 

1.10 

0.8313 

0.90 

6 

1193.9 

0.0130 

5.18 

2.20 

0.6004 

1.00 

7 

1193.5 

0.0261 

5.00 

4.39 

0.1922 

1.10 

8 

1192.6 

0.0533 

4.72 

8.18 

0.0243 

1.33 

9 

1196.0 

0.0049 

6.35 

1.10 

0.9916 

1.00 

10 

1195.9 

0.0097 

6.25 

1.96 

0.7091 

0.97 

11 

1195.9 

0.0195 

6.26 

4.03 

0.3293 

1.03 

12 

1196.1 

0.0395 

6.45 

7.81 

0.0727 

1.10 

16 

1193.9 

0.0065 

5.15 

0.97 

0. 9442 

0.84 

17 

1193.8 

0.0130 

5.10 

2.20 

0.5683 

1.03 

20 

597.1 

0.0130 

5.28 

2.32 

0.8255 

0.96 

21 

596.9 

0.0261 

5.15 

4.39 

0.5759 
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APPENDIX  B 

HORIZONTAL- VERTICAL  FORCE  INTERACTION  AT  SMALL  HORIZONTAL 

INCLINATION  ANGLES 


Equations  (37)  and  (38)  may  be  combined  to  give 


(F„) 

N  x 


2  2 
cos  4> ( 1  -  cos  q  cos 

"  sin  0 


cos  9  sin^  0 


(B-I) 


A  typical  operational  scenario  would  have  9  vary  as  a  zero  mean  random 
variable  in  both  space  and  time  with  a  maximum  value  of  less  than  3 
degrees.  The  depression  angle  <b  would  have  a  steady  value  between  5  and  45 
degrees.  For  this  range  of  4  and  4>  we  may  approximate  equation  (B-l)  as 


f  )  *  -R  |~- 

N  x  j 


(B-2) 


The  interpretation  in  the  text  is  that  the  bracketed  term  of  equation 
(B-2)  serves  to  modify  (amplify)  R^  or,  more  appropriately,  Cjj  from 
equation  (34). 

The  point  is  that  the  singular  behavior  of  the  amplifying  term  near  9 
equal  to  zero  does  not  imply  a  singular  behavior  in  (Fjj)x,  but  rather  a 
modification  of  its  horizontal  inclination  angle  dependence  from  a  92  to  a 
9  dependence.  The  singularity  of  equation  (38)  is  the  consequence  of  some 
mathematical  "bookkeeping"  resulting  from  that  functional  change. 

To  avoid  the  singularity  in  question  formally,  we  could  define  an 
effective  amplification  (Agpp)  as 


V,] 

Vxl<j>-ol 


(B-3) 


where  E  is  the  expectation  operator  and  9  is  a  zero  mean,  normally 
distributed  random  variable  with  a  root  mean  square  value  of  aQ.  The 
numerator  and  denominator  of  equation  (B-3)  are  evaluated  from  the  integrals 


E  [(v,U-o] 


-2R_  r  2e-l/2(8/0<;)‘  „„ 

2tt  a  J  o 
o 


(B-4) 


[( V.  ] 


2 

-2R  cos  sin 


,.-i/2(e/o0)‘  . 


\  ■  ?".•  V:  V 
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The  integrals  of  equation  (B-4)  may  be  found  in  integral  tables.  The 
results  of  the  integration  are  then  substituted  into  equation  ( B— 3 )  to  give 


s^  $  sin  <j> 


(B-5) 


A  comparison  of  equations  (38)  and  (B-5)  shows  them  to  be  functionally  the 
same;  yet  with  (B-5)  we  have  avoided  the  singular  problem. 
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